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The small-x contributions to the Bjorken sum rule within unified picture 
ln 2 x+LO DGLAP for different input parametrisations g^ s (x, Qq) are pre- 
sented. Theoretical predictions for J ' gi S (x, Q 2 = 10)dx are compared 
with the SMC small- a; data. Rough estimation of the slope A, controlling 
the small- a: behaviour of gf s ~ x~ x from the obtained results and SMC 
data is performed. The crucial role of the running coupling a s = a s (Q 2 / z) 
at low-x is taken into account. 

PACS numbers: 12.38 Bx 

1. Introduction 

The results of SIDIS (semi inclusive deep inelastic scattering) experi- 
ments with polarised beams and targets enable the extraction of the spin 
dependent quark and gluon densities. This powerful tool of studying the 
internal spin structure of the nucleon allows verification of sum rules. One 
of them is the Bjorken sum rule (BSR) p], which refers to the first moment 
of the nonsinglet spin dependent structure function gi S (x, Q 2 ). Because of 
SUf(2) flavour symmetry, BSR is regarded as exact. Thus all of estimations 
of polarised parton distributions should be performed under the assumption 
that the BSR is valid. Determination of the sum rules requires knowledge 
of spin dependent structure functions over the entire region of x £ (0; 1). 
The experimentally accessible x range for the spin dependent DIS is how- 
ever limited (0.7 > x > 0.003 for SMC data [2]) and therefore one should 
extrapolate results to x = and x = 1. The extrapolation to x — > 0, where 
structure functions grow strongly, is much more important than the extrap- 
olation to x — ► 1, where structure functions vanish. Assuming that the 
BSR is valid, one can determinate from existing experimental data the very 
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small- x contribution (0.003 > x > 0) to the sum rule. Theoretical analysis 
of the small-x behaviour of gi S (x,Q 2 ) = g\{x,Q 2 ) — g™(x,Q 2 ) together 
with the broad x-range measurement data allow verification of the shape 
of the input parton distributions. In this way one can determinate the free 
parameters in these input distributions. Experimental data confirm the 
theoretical predictions of the singular small-x behaviour of the polarised 
structure functions. It is well known, that the low-x behaviour of both 
unpolarised and polarised structure functions is controlled by the double 
logarithmic terms (a s ln 2 x) n For the unpolarised case, this singular 

PQCD behaviour is however overridden by the leading Regge contribution 
[3]. Therefore, the double logarithmic approximation is very important par- 
ticularly for the spin dependent structure function g\ . The resummation 
of the ln 2 x terms at low x goes beyond the standard LO and NLO PQCD 
evolution of the parton densities. The nonsinglet polarised structure func- 
tion gi S , governed by leading a™ln 2n x terms, is a convenient function both 
for theoretical analysis (because of its simplicity) and for the experimental 
BSR tests. The small-x behaviour of g± s implied by double logarithmic 
approximation has a form x~~ A with A ~ 0.4. This or similar small-x ex- 
trapolation of the spin dependent quark distributions have been assumed in 
recent input parametrisations e.g. in [5] . [7| . [T3~] . More singular parametri- 
sation of gi NS (x,Ql = 4) ~ x~ a8 at small-x, based on the QCD (LO 
and NLO) analysis of the world data on polarised deep inelastic scattering, 
has been presented in |14j . Mentioned above double logarithmic approach 
is however inaccurate for QCD analysis at medium and large values of x. 
Therefore the double logarithmic approximation should be completed by 
LO DGLAP Q 2 evolution. In our theoretical analysis within /n 2 x+LO 
DGLAP approach we estimate g^ s at low-x and hence the small-x con- 
tributions f*° g? s (x,Q 2 )dx, \ll 9i S (x,Q 2 )dx (x ,xi,x 2 <. 1) to the BSR 
for different input quark parametrisations: the Regge nonsingular one and 
the singular one. We compare our results with the suitable experimental 
SMC data for BSR. In the next section we recall some of the recent the- 
oretical developments concerning the small-x behaviour of the nonsinglet 
polarised structure function gi S . Section 3 is devoted to the presentation 
of the unified ln 2 x+LO DGLAP approximation. We also discuss the role 
of the running coupling a s . Section 4 contains our results for the struc- 
ture function g± s at small-x and for contributions to the Bjorken sum rule 
AIbsr(xi,X2,Q 2 ) = III 9i S (x,Q 2 )dx (xi,x 2 < 1). We present our pre- 
dictions using flat (nonsingular) ~ (1 — x) 3 and singular ~ x~ x at small-x 
parametrisations of the input structure function gf^x,^) as well. We 
compare our results with the SMC data for the small-x contribution to the 
BSR. We roughly estimate the slope A controlling the small-x behaviour 
of g^ s ~ x~ x from our g^ s predictions and from the SMC data, basing 
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on the validity of the BSR. We compare also results AIbsr(x\,X2, Q 2 ) and 
g^ s (x = 10~ 6 ,(3 2 = 10) in different approximations: pure LO DGLAP, 
pure ln 2 x, In 2 x+LO DGLAP and obtained for different a s parametrisa- 
tions: a s = const, a s = a s (Q 2 ), a s = a s (Q 2 /z). Finally, Section 5 contains 
a summary of our paper. 

2. Small- a; behaviour of the nonsinglet spin dependent structure 

function g™ s (x,Q 2 ) 

The small value of the Bjorken parameter x, specifying the longitudinal 
momentum fraction of a hadron carried by a parton, corresponds by defini- 
tion to the Regge limit (x — > 0). Therefore the small- x behaviour of struc- 
ture functions can be described using the Regge pole exchange model (5). In 
this model the spin dependent nonsinglet structure function g^ s = g\ — Q\ 
in the low-x region behave as: 

g? s (x,Q 2 )=~f(Q 2 )x-^ (2 .i) 

where 0^(0) is the intercept of the A\ Regge pole trajectory corresponding 
to the axial vector meson and lies in the limits 

- 0.5 < ctAi (°) < (2.2) 

This low value of the intercept (|2.2[) implies the nonsingular, flat behaviour 
of the g^ s function at small-x. The nonperturbative contribution of the A\ 
Regge pole is however overridden by the perturbative QCD contributions, 
particularly by resummation of double logarithmic terms ln 2 x. In this way 
the Regge behaviour of the spin dependent structure functions is unstable 
against the perturbative QCD expectations, which at low-x generate more 
singular x dependence than that implied by l|2.1|) - (|2.2|) . Nowadays it is 
well known that the small- x behaviour of the nonsinglet polarised structure 
function g± s is governed by the double logarithmic terms i.e. (a s ln 2 x) n 
0,01 . Effects of these ln 2 x approach go beyond the standard LO and 
even NLO Q 2 evolution of the spin dependent parton distributions and 
significantly modify the Regge pole model expectations for the structure 
functions. From the recent theoretical analyses of the low-x behaviour of 
the gi S function |5] one can find that resummation of the double logarithmic 
terms (a s ln 2 x) n leads to the singular form: 

g? s (x,Q 2 )~x- x (2.3) 

with A ~ 0.4. This behaviour of g? s is well confirmed by experimental 
data, after a low-x extrapolation beyond the measured region [2] . [TU] . [TT] . 
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3. Unintegrated structure function f NS (x,Q 2 ) within double 
logarithmic ln 2 x and unified ln 2 x-\-hO DGLAP 
approximations 

Perturbative QCD predicts a strong increase of the structure function 
g^ s (x, Q 2 ) with the decreasing parameter x 0)11] what is confirmed by 
experimental data P].|1U|.[TT]. This growth is implied by resummation of 
ln 2 x terms in the perturbative expansion. The double logarithmic effects 
come from the ladder diagram with quark and gluon exchanges along the 
chain. In this approximation the unintegrated nonsinglet structure function 
/ (x, Q 2 ) satisfies the following integral evolution equation 0: 

i Q 2 /z /2 

f NS (x,Q 2 ) = fo NS (*) + [- I ^a s f NS ^,k' 2 ) (3-1) 



where 



z J k' 2 

x o2 



* = ^ (3.2) 
and /q (x) is a nonperturbative contribution which has a form: 

f» s (x) = a s }%r S (z) (3-3) 

gi (x) is an input parametrisation 

9° 1 NS W=9? S (x,Q 2 = Q 2 o) (3-4) 
The unintegrated distribution f NS (x,Q 2 ) is related to the g\ (x,Q 2 ) via 

Eq. p.lj) generates the leading small-x behaviour of f NS and hence g^ s , 

but it is inaccurate in describing the total Q 2 evolution. For larger values 

l 

of x, which are involved in the evolution equation ()3.1|) via J dz one should 

X 

take into account Q 2 DGLAP evolution with complete splitting function 
Pq q (z). Therefore the double logarithmic approach should be completed 
by LO DGLAP Q 2 evolution. Unified description of the polarised struc- 
ture function f NS (x,Q 2 ) incorporating DGLAP evolution and the double 
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logarithmic ln 2 x effects at low- a; leads to the following equation for the 
unintegrated distribution f NS (x,Q 2 ) [T8| : 



l Q 2 /z 

ii.'- 



f NS (x,Q 2 ) = f NS (x) + J *1 J ^a s f^,k 



dk' 2 fdz (1 + z 2 )f(x/z,k' 2 ) -2zf(x,k' 2 ) 
+ I -nr- / — a s 



k' 2 J z 1 — z 

o 2 

J dk' 2 ,3 



f Hie'* 3 

+ *.J ^(| + 21n(l -*))/(*, A;' 2 ) (3.6) 



where 



) dz{l + z 2 )gf\x/z)-2zgf\x) 



z 1 — z 



x 



+ (l + 2ln(l-x))g^(x)] (3.7) 

The unintegrated distribution f NS in the equation (|3.6|) is related to the 
< s (z,Q 2 )via 

Q 2 (l/x-l) 

r dk 2 / k 2 ^ 



An important role in solutions of Q3.1|) and (|3.6|) plays the coupling a s , 
which can be parametrised in different way. The simplest choice of a s is a 
Constance (nonrunning) coupling: 

a s = const (3-9) 

This simplification allows the analytical analysis of the suitable evolution 
equations for truncated and full moments of the unintegrated structure func- 
tion f NS (x,Q 2 ) within ln 2 x approximation [H],[E]- The introduction of the 
running coupling effects implies a s in ()3.1|) and ()3.6|) of a form 

a s = a s (Q 2 ) (3.10) 

It has been however lately proved [S], that dealing with a very small- x 
region one should use a prescription for the running coupling in a form a s = 
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a s (Q 2 /z). This parametrisation is theoretically more justified than a s = 
a s (Q 2 ). Namely, the substitution a s = a s (Q 2 ) is valid only for hard QCD 
processes, when x ~ 1. However the evolution of DIS structure functions at 
small- 2; needs "more running" a s : 

a s = a s (Q 2 /z) (3.11) 

Our predictions for g^ s and AIbrs(xi, £2, Q 2 ) for different forms of a s will 
be presented in the forthcoming section. 



4. Predictions for g^ s and small- 2: contribution to the BSR 

Our purpose is to calculate the nonsinglet polarised structure function 
g^ s (x, Q 2 ) and hence also the contribution to the Bjorken sum rule in the 
small- x region. The BSR is a fundamental rule and must be hold as a 
rigorous prediction of QCD in the limit of the infinite momentum transfer 
Q 2 : 

1 

W = r?-r?= ( dxg? s (x,Q 2 ) = \\ 9 -^\ (4.1) 
J 6 g v 



where 



r? = y dxgl{x,Q 2 ) (4.2) 





1 



1 1 



= I dxg™(x,Q 2 ) (4.3) 



and |^| is the neutron /3-decay constant 

|M| = F + D = 1.2670 (4.4) 
9V 

Hence the BSR for the flavour symmetric sea quarks scenario (Au = Ad) 
reads: 

1 

Ibsr(Q 2 ) = J dxgf s (x, Q 2 ) w 0.211 (4.5) 


The small-x contribution to the BSR has a form: 

AI BSR (x 1 ,x 2 ,Q 2 ) = Jdxg? s (x,Q 2 ) (4.6) 
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Below we present our results for g^ s and AIbsr at small- a; obtained for 
different a s sets (j3.9p - p.lip within combined ln 2 x+LO DGLAP approach. 
We compare these predictions with pure LO DGLAP and pure ln 2 x results 
as well. We solve numerically the evolution equation ()3.6p in a case of unified 
ln 2 x+LO DGLAP picture and in a case of pure LO DGLAP, when one gets 
the following equation: 

f NS (x,Q 2 ) = f NS (x) + 

Q dk' 2 } dz {1 + z 2 )f(x/z,k' 2 ) -2zf(x,k' 2 ) 
-a* 



k' 2 J z 1 — z 

dk' 2 ,3 



f dk' i 3 



In order to have comparable results, for pure ln 2 x analysis we also use nu- 
merical solutions of ()3.ip . Our predictions have been found for two different 
input parametrisations g® (a;), chosen at Qq = lGeV 2 : 



1. fli (x) = 0.8447(1 - xy (4.8) 

2. g^ NS (x) = 0.290x- a4 (l-x) 25 (4.9) 

Input 1 is the simple Regge form, Constance as x — > 0; input 2 is a "toy" 
model, in which we have used the latest theoretical results concerning the 
small- x behaviour x -0,4 of the nonsinglet function gi S [SJ. In Fig.l we plot 
inputs g1 NS (x) 1)4.811 - l)4.9p in the low-x region [10~ 5 10 -2 ] together with 
g^ s (x,Q 2 ) results for Q 2 = lOGeV 2 within /n 2 x+LO DGLAP approach 
with "very running" coupling a s = a s (Q 2 /z). Fig. 2 contains comparison of 
three approximations: pure ln 2 x, pure LO DGLAP and unified ln 2 x+LO 
DGLAP. We present structure function g^ s (x, Q 2 = 10) for a s = a s {Q 2 / z). 
Finally, Fig. 3 shows the g^ s (x, Q 2 = 10) results within combined Zn 2 x+LO 
DGLAP approach for different parametrisations of a s : (|3.9p - (|3.1ip . In each 
figure we present the solutions for both input parametrisations (|4.8p - (|4.9|) . 
Numbers at each plot correspond to the suitable inputs 1 or 2. In Table I 
we present our results for the low-x contributions to the BSR (|4.6p together 
with e{x\,X2), which is defined by the following expression: 

J dxg? s (x, Q 2 ) = [1 + e{x 1 ,x 2 )} J dxgl NS {x) (4.10) 
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Fig. 1. Input parametrisations g1 NS (dashed) and g^ s (x,Q 2 = 10) 

(solid) for these inputs within ln 2 x+LO DGLAP approach and for running 
a s (Q 2 /z). 
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Fig. 2. The small- a; predictions for g^ s (x,Q 2 = 10) within different approxima- 
tions: LO DGLAP (dotted), ln 2 x (dashed), unified ln 2 x+LO DGLAP (solid). 
Plots for both inputs l|4.8|) - (|4.9|) and running a s (Q 2 /z). 



In the last column we give the percentage value p[ 



P = AlB f {xi ;^; Q2) -m% (4.ii) 

Ibsr{Q z ) 
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Fig. 3. The small- a; predictions for gi S (x,Q 2 = 10) within unified /n 2 a;+LO 
DGLAP approach for different a s : a s = 0.18 (dotted), a s {Q 2 ) (dashed), a s (Q 2 / ' z) 
(solid). Plots for both inputs (gj^-fOj!. 



In Table II we collect results for all possible combinations of approximations 
and a s sets: LO DGLAP and a s = const = 0.18, LO DGLAP and a s = 
a s (Q 2 ),..., ln 2 x and a s = a s (Q 2 / z) etc. We present here AIbsr(0, 3 • 
10 _3 ,10), gi s (x = W~ 6 ,Q 2 = 10). In the last column the effective slope 
A at Q 2 = lOGeV 2 and small-x [10- 6 ;10- 5 ] is shown. We use again 

both inputs (x). 



Xl 


X2 


AI B sr(xi,x 2 , 10) 


e(xi,x 2 ) 


p% 





3 • 10~ 3 


(1) 0.006108 


1.4213 


2.89 


(2) 0.020668 


0.4097 


9.80 





10- 2 


(1) 0.016050 


0.9289 


7.61 


(2) 0.040000 


0.3336 


18.96 


10" 5 


10" 3 


(1) 0.002457 


1.9422 


1.16 


(2) 0.010450 


0.4574 


4.95 


10" 4 


io- 2 


(1) 0.015672 


0.9028 


7.43 


(2) 0.037380 


0.3214 


17.72 



Table 1. The small-x contribution to the BSR (|4.6(l for different input parametrisa- 
tions (|4.8|) - (|4.9|) within ln 2 x+LO DGLAP approximation with running a s (Q 2 /z). 

From these results one can read that the low-x g^ s (x, Q 2 ) values and hence 
the low- a; contributions to the BSR strongly depends on the input parametri- 
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approach 


a s 


A/ B5/? (0, 0.003) 


N 4 
01 


A 


l.Regge 


LO 


const=0.18 


0.003879 


2.07 


0.06 


a s (Q z ) 


0.005742 


4.31 


0.11 


a s (Q z /z) 


0.004534 


2.17 


0.04 


ln 2 x 


const=0.18 


0.005614 


10.2 


0.25 


a s (Q z ) 


0.008855 


23.9 


0.31 


a s (Q z /z) 


0.007043 


9.93 


0.20 


ln 2 x+LO 


const=0.18 


0.005440 


9.75 


0.25 


a s (Q' z ) 


0.008281 


21.8 


0.30 


a s (Q z /z) 


0.006108 


7.36 


0.18 


2." Toy" 


LO 


const=0.18 


0.017283 


88.6 


0.40 


a s (Q z ) 


0.020543 


110 


0.40 


a s (Q z /z) 


0.019184 


98.0 


0.40 


ln 2 x 


const=0.18 


0.019026 


113 


0.42 


a s (Q z ) 


0.023641 


161 


0.44 


a s (Q z /z) 


0.022175 


130 


0.41 


ln 2 x+LO 


const=0.18 


0.018756 


111 


0.42 


a s (Q z ) 


0.022712 


152 


0.43 


a s (Q z /z) 


0.020668 


117 


0.41 



Table 2. AI BSR (0, 0.003, 10), g^ s (x = 10" 6 ,10) and A for both input parametri- 
sations (|4.8|l - l|4.9(l within different approaches and a s . 



sation g^ NS . For the flat Regge form (|OJ) AI BS r{0, 10~ 2 , 10) is equal 
to around 7.6% of the total Ibsr = 0.211, while for the singular input 
(l43|) 19.0%. The structure function #f s itself at very small-x = 10 6 and 
Q 2 = lOGeV 2 is in a case of x~ 0A input about 16 times larger than for the 
flat one. The value of e(xi,X2), defined in (|4.10j) varies from 0.9 -r- 1.9 for 
the Regge input 1 to 0.3 4- 0.5 for the singular input 2. The effective slope A 
(|2.3|) describing the small- x behaviour of the structure function g^ s remains 
unchanged in a case of the singular input. Namely, the x~ 0A shape of the 
input Qi s (x) implies again the same low-j; behaviour of the g^^x, Q 2 ), 
independently of the Q 2 -evolution approach. Quite different situation oc- 
curs for the flat inputs e.g. the Regge one (|4.8|) . where the singular small-x 
behaviour of the gi S (x, Q 2 ) is totally generated by the QCD evolution with 
ln 2 x terms. Pure double logarithmic ln 2 x approach or combined ln 2 x+LO 
DGLAP approximation give the value of A from 0.2 to 0.3. Only in the 
pure LO DGLAP analysis we obtain A < 0.1. It means that the double 
logarithmic ln 2 x effects are better visible in a case of nonsingular inputs. 
In a case of singular input parametrisations g\ NS ~ x~ x (e.g. A ~ 0.4) the 
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growth of g^ s at small-x, implied by the ln 2 x terms resummation, is hidden 
behind the singular behaviour of g® NS , which survives the QCD evolution. 
Comparing plots from Fig. 2 and results from the last column in Table II, 
one can read that the ln 2 x resummation gives steep growth of the g^ s in 
the small-x region. It is well visible in a case of the Regge input, where 
gr^ (x, Q 2 ) within ln 2 x or /n 2 x+LO DGLAP approaches strongly dominate 
over that, obtained in pure LO DGLAP approximation. Double logarith- 
mic contributions of the type (a s ln 2 x) n , which lead to the strong growth 
of structure functions at low-x are not included in the DGLAP evolution 
(LO or NLO). Differences between pure ln 2 x and In 2 x+LO DGLAP results 
within the same set of a s are not very significant. However, pure ln 2 x ap- 
proximation overestimate the value of gi S and should be accompanied by 
the LO DGLAP evolution. This is because in the larger- x region, involved 
in the evolution equation for f NS (|3.1|) . pure ln 2 x analysis is inadequate. 
The crucial point in QCD analysis is a treatment of the coupling a s . The 
problem of a s parametrisations in high energy processes has been widely 
discussed in [S]. Fixed, Constance a s is very convenient in many physical 
problems. Thus, use of the fixed a s simplifies the evolution equations for 
structure functions and enables easy analytical solutions. However, an in- 
troduction of the fixed coupling needs a reasonable scale and this is somehow 
"artificial". Namely, the scale for fixing of a s is not well defined. In pertur- 
bative QCD one should take into account running a s effects. In this way, 
usually, the prescription for the running coupling reads a s = a s (Q 2 ). This 
construction is however accurate only for hard QCD processes, where x ~ 1. 
On the other hand, many interesting QCD processes (e.g. DIS at low-x) are 
Regge-like. For these cases, with small-x involved, "hard" running a s (Q 2 ) 
is incorrect. Instead of a s (Q 2 ) one has to use a modified parametrisation of 
a s : 

a s = a s (k 2 J(3) (4.12) 

where k\ is the transverse momentum of the ladder parton and (3 is the 
standard Sudakov parameter. In our approach this prescription reads as 
(|3.11j) . From Fig. 3 and Table II we are able to compare the predictions 
for g± s at small-x for three a s parametrisations (|3.9|) - (|3.11j) . The results 
for a s = const = 0.18 and a s = a s (Q 2 /z) (within the same approach LO 
or ln 2 x etc. and with the same input g® NS ) are similar but significantly 
smaller than in a case of running a s = a s (Q 2 ). Within In 2 x+LO DGLAP 
approximation with flat input for "very" running a s (Q 2 /z) (x < z < 1), via 
weaker coupling {a s (Q 2 / z) < a s (Q 2 )), the value of gi S (x = 1CT 6 , Q 2 = 10) 
is almost 3 times smaller than for the "hard" running a s (Q 2 ). It is a good 
lesson how choice of the running coupling influences the results in the low-x 
region. From the experimental SMC data [TO] the low-x contribution to the 
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BSR at Q 2 = lOGeV 2 is equal to 



0.003 



6 J g? s (x,Q 2 = 10)dx = 0.09 ±0.09 (4.13) 
o 

The above result has been obtained via an extrapolation of g^ s to the un- 
measured region of x: x — > 0. Forms of the polarised quark distributions 
have been fitted to SMC semi-inclusive and inclusive asymmetries. In the 
fitting different parametrisations of the polarised quark distributions |15j 
|16j have been used. The extrapolation of g^ s to very small-x region de- 
pends strongly on the assumption (input parametrisation) made for this 
extrapolation. In this way present experimental data give only indirectly 
the estimation of the small-x contribution to the moments of parton dis- 
tributions. The result (|4.13|) with a large statistical error and strongly 
fit-dependent cannot be a final, crucial value. Nevertheless we would like 
to estimate the exponent A in the low-x behaviour of gi S ~ x~ A using the 
above SMC result for the small-x contribution to the BSR. Assuming the 
validity of the BSR (gSJ) at large Q 2 = lOGeV 2 , one can find: 

x 1 

J dxg? s (x, Q 2 ) = Ibsr{Q 2 ) ~ J dxg? s (x, Q 2 ) (4.14) 

x 

where xo is a very small value of the Bjorken variable. Taking into ac- 
count the small-x dependence of gi S ~ x~ x and the experimental data for 
AIbsr(0, 0.003, 10) one can obtain: 

0.003 

C J x~ x dx = 0.015 ±0.015 (4.15) 
o 

The constant C can be eliminated from a low-x SMC data jlOj : 

Cx~ x = g™~ p (x, 10) (4.16) 

Taking different small-x SMC data, we have found A = 0.37 (x = 0.014); 
A = 0.20 (x = 0.008); A = 0.38 (x = 0.005). 

It seems nowadays that the most probably small-x behaviour of g± s is 

5 f s (x,Q 2 )~x- ' 4 (4.17) 



This results from latest theoretical analyses [5] , which take into account the 
running coupling effects at low-x. It has been shown in [S] that the intercept 
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A controlling the power-like small-x behaviour of g± depends on the choice 
of the parameters rif (flavour number), Qq (input scale) and Aqcd- The 
maximal value of A, which gives the maximal contribution to the structure 
function g? s in the perturbative QCD, incorporating ln 2 x effects, is equal to 
0.4. The same value of A = 0.4 was obtained in the semi-phenomenological 
estimation from BSR for lower Q 2 |17) . The polarised (nonsinglet and singlet 
as well) structure functions are presently the objects of intensive theoretical 
investigations. Maybe the crucial point for understanding of the small-x 
behaviour of structure functions is an analysis beyond the leading order 
a™ln 2n x. The resummation of a^ n+1 ' ) /n 2n x terms is studied in j^Uj. The 
corrections to gi S due to the nonleading terms are on the level of 1% in 
the accessible at present experimentally x region, but can be larger (up to 
about 15%) at very small-x ~ 10~ 5 . In the situation, when the present 
experimental data do not cover the whole region of x G (0; 1), theoretical 
predictions for e.g. structure functions in the unmeasured low- a; region 
cannot be directly verified. Latest experimental SMC 0, ^U] and HERMES 
jllj data provide results for the BSR from the region 0.003 < x < 0.7 and 
0.023 < x < 0.6 respectively. In the very small-x region exist only indirect, 
extrapolated results with large uncertainties. Small-x contribution to the 
Bjorken sum rule resulting from such indirect SMC data analysis is equal to 
0.015±0.015. Large uncertainties of the small-x experimental results disable 
unfortunately realistic comparison the data with the theoretical predictions. 
Namely, all our results for AI BS r(0, 0.003, 10) in Table II: from 0.004 (for 
LO, a s = 0.18, flat input) to 0.024 (for ln 2 x, a s (Q 2 ), singular x~ 0A input) 
are in agreement with SMC data (within the total error). Nevertheless, the 
progress in theoretical 0,01, , [IE] , d , |H1 and experimental [2]. [Tn].[TT] 
investigations gives hope that our knowledge about structure functions at 
small-x is getting better. 



5. Summary and conclusions 

In this paper we have estimated the nonsinglet polarized structure func- 
tion g^ s at small-x and also contributions from the small-x region to the 
Bjorken sum rule. We have used the numerical solutions within unified 
double logarithmic and DGLAP (ln 2 x+LO DGLAP) approximation. Our 
predictions for g± (x,Q 2 ) and AIbsr{xi, %2, Q 2 ) have been found for two 
input parametrisations g1 NS (x, Qq). These parametrisations describe dif- 
ferent small-x behaviour of g® NS = at Qq: g^ s ~ x~ x . The main 
conclusion from our analyses is that the structure function g^ s at small-x 
and hence also the small-x contribution to the BSR strongly depends on the 
input parametrisation gi NS . The percentage value AIbsr{0, 10 -2 , Q 2 = 10) 
of the total BSR ~ 0.211 varies from 7.6 for the flat Regge input 1 (A = 0) to 
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almost 19 for the singular one 2 (A = 0.4). The structure function g^ itself 
at very small-x = 1CT 6 and Q 2 = lOGeV 2 is in a case of x~ 0A input about 
16 times larger than for the flat one. Double logarithmic ln 2 x effects, re- 
sponsible for the strong growth of the structure function in the low-x region, 
are better visible in a case of nonsingular inputs. In a case of singular input 
parametrisations g® NS ~ x~ A (e.g. A ~ 0.4) the growth of g^ s at small- 
x, implied by the ln 2 x terms resummation, is hidden behind the singular 
behaviour of g® NS , which survives the QCD evolution. Input parametrisa- 
tion 2 incorporates latest theoretical investigations, which suggest singular 
small-x shape of polarised structure functions: ~ x~ 0,4 for the nonsinglet 
case and even ~ x -0 - 8 for the singlet one. Both these values are indirectly 
confirmed by fitted experimental HERMES data. Basing on these results, 
similar extrapolations of the spin dependent quark distributions towards the 
very low-x region have been assumed in several recent input parametrisa- 
tions Aq(x,Qo). Our results for the small-x contribution < x < 0.003 to 
the BSR are in agreement with the experimental SMC data (for both in- 
puts). However it must be emphasized, that SMC data for the low-x region 
suffer from large uncertainties. Using SMC data for gi S at small-x (0.14, 
5 • 10~ 3 , 8 • 10~ 3 ) we have estimated the exponent A which governs the low-x 
behaviour of g^ s . Thus we have obtained A = 0.20 -j- 0.38 with large uncer- 
tainties. This effective slope A calculated for low-x € [10~ 6 ; 10~ 5 ] from g^ s 
in our approach amounts about 0.2 (for running a s (Q 2 / z) and Regge-like 
flat input). In order to have reliable theoretical predictions for the polarised 
structure function g^ s (x, Q 2 ) we have used unified approach which contains 
the resummation of the ln 2 x and the LO DGLAP Q 2 evolution as well. It is 
because the pure ln 2 x approximation generates correctly the leading small- 
x behaviour of the polarised structure function but is inaccurate for larger 
values of x. Another crucial point of the presented analysis is the role of the 
running coupling effects. Latest theoretical studies suggest introduction of 
the running coupling of a form a s = a s (Q 2 / z) instead of a s = a s (Q 2 ). This 
is more justified in the small-x region. We have found that the choice of 
the running coupling significantly influences the results in the low-x region. 
E.g. the value of g? s (x = 10~ 6 ,Q 2 = 10) is for "hard" running a s {Q 2 ) 
almost 3 times greater than for the "very" running a s {Q 2 / z) (in a case of 
nonsingular input gi NS ). Proper theoretical treatment of the Q 2 evolution 
of structure functions in the whole (small and large) x region with all essen- 
tial perturbative leading and even nonleading effects involved should be a 
subject of further intensive investigations. It is important because of lack of 
the experimental data from the very small-x region (x < 0.003). Agreement 
of the theoretical predictions e.g. for the BSR with real experimental data 
at medium and large x may give hope, that for the very interesting small-x 
region the suitable theoretical results are also reliable. 
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